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ABSTRACT 


This study was conducted at the Naval Postgraduate School 
to investigate the formation of particle flow cells at mini- 
mum fluidization gas flow rates in a rectangular gas-fluidized 
bed. The primary objective is to determine if multiple cell 
formation occurs as the bed width to height ratio increases. 
The settled bed height iS Maintained at a constant level for 
all runs. The secondary objective is to determine if the 
presence of external plate heaters at the end walls of the 
particle bed will alter the particle flow cell formation 
patterns when energized and to what extent they affect the 
Pearae mixing (heat transfer) flow patterns within the par- 
ticle flow. Pressure measurements across the bed are used 
to determine minimum fluidization points and thermocouple 
probe measurements determined the thermal gradients within 
the bed for each bed configuration. Visual observations on 
the external bed movement are included in the determination 


of the results of this experimental study. 
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NOMENCLATURE 


English Letter Symbols 


a Constant, dependent on tube orientation - Reference 30. 

A. Cross-sectional area of bed container. 

A. Contact area of packet and surface - Figure 18. 

B Shape factor - Reference l2. 

Bi Biot number. 

Ca Specific heat of packet or quiescent bed. 

Sot Spectiite Neat of fluid. 

Cos Seecitic neat of solid particle. 

Cp Constant, correction for location of immersed body - 
Reference 34. 

C. Speecliic Meat Of Solid particle. 

dp Diameter of sphere of same specific area as that 
of the bed particle. 

D Diameter of the fluidized bed. 

D- Effective diameter of free area across the bed. 

Dis Diameter of an immersed heater. 

D. Particle diameter of sieved fraction. 

5 Particle diameter. 

Dy Average particle diameter. 

D. Distance between two successive layers - Figure 16. 

Di Diameter of immersed tube. 

D5 Diameter of immersed tube being equal to 20 mm. 

g Acceleration of gravity. 


iE 





Ca 


Ga 


Fluid mass velocity. 
Fluid mass velocity based on voidage area. 
Galileo number (Archimedes number). 


Hyeoseiecalcal ninva Mass velocity required to 
merely expand a bed of particles. 


Fluid mass velocity through a fluidizing bed. 
Fluid mass velocity for minimum fluidization. 


Optimum fluid mass velocity giving a maximum 
value of the heat transfer coefficient. 


Fluid mass velocity based on kinematic viscosity. 


Heat transfer coefficient for the contact 
resistance at the wall. 


Empirical constant - Reference 44. 


Heat transfer coefficient between the fluidized 
bed and the surface. 


Instantaneous local heat transfer coefficient. 
Local average heat transfer coefficient. 


Maximum heat transfer coefficient between 
fluidized bed and surface. 


Superficial bubble-to-cloud heat transfer coefficient. 
uid thermal conductivity . 

Thermal conductivity of packet or quiescent bed. 
maermnal Conductivity of particle. 

Thermal conductivity of sphere. 

Distance from bed top in packet theory - Reference 47. 
Fluidized bed height. 
Expanded bed height. 


Length of heat exchange surface. 


dg 





Beaenotonenat DOlne Of Minimum fluidization. 
SlOpeaeeeigure 12. 

Slope - Figure lL5. 

Fluidization number. 

Nusselt number. 

Maximum value of Nusselt number. 
Empirical constant - Reference 44. 
Center-to-center distance (pitch) of adjacent tubes. 
Prandtl number. 


Heat flux. 


Heat transferred into the core portion of the 
an interchange of particles - Figure 23. 


bed by 


Heat flow rate from the wall into a packet. 
Instantaneous heat transfer rate. 

Heat taken by particles flowing parallel to the 
SsuUbEace in a second zone Of emulsion of thickness 


5 er heise 2 0 


Fluctuation ratio - highest to lowest bed level ina 
gas-fluidized charge. 


Radius of sphere. 

Bed expansion ratio. 
Reynolds number. 
Reynolds number based on particle free fall velocity. 
Modified Reynolds number - Reference 27. 
Reynolds number at minimum fluidization. 


Optimum value of Reynolds number giving a maximum 
heat transfer coefficient. 


Reynolds number based on bed diameter. 


aS 





Greek 


a 


Average replacement of packets at wall by side mixing. 
Stanton number. 

Time. 

Bulk bed constant temperature. 

Fluid film temperature. 

Gas temperature into bed. 

Gas temperature out of bed. 

Temperature of wall or heater. 

Pemeraleiiitd velocity with no particle motion. 
Fluid velocity required for initial expansion. 
Fluid velocity required for minimum fluidization. 
SMpecrenciakn velocity at tluidization conditions. 
Particle velocity. 

Pita velocity during fluidization. 

Particle velocity of Packet Model. 

Merminal fluid velocity. 

Volume displaced by particles. 

Bed width. 

Axis direction - Figures 21 and 22. 

Weight fraction of sieved particles for D.. 
Packet width - Figures 21 and 22. 

Horizontal pitch of tube bundle - Reference 47. 


Vertical pitch of tube bundle - Reference 47. 


Letter Symbols 


Constant - Reference 9@. 
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Constant - Reference 9. 
Specific weight of fluid. 


Boundary layer thickness at point between two particle 
layers - Figure 16. 


Average boundary layer thickness at point between two 
points - Figure 16. 


Thickness of emulsion layer. 
Thickness of gas layer. 


Minimum distance between surfaces of adjacent tubes 
(gap) . 


Pressure drop across bed at Ug: 
Pressure drop across bed at U.: 
Pressure drop across bed at Une? 


Ideal pressure drop across the bed calculated from the 
buoyant weight of the bed. 


VOlGurtacrei1on at fluidization conditions. 
Initial particle voidage as filled. 
Particle voidage during expansion. 
Particle voidage at minimum fluidization. 
Poleamitaction at fluidization conditions. 
Eeutarzaction CLEICIency . 

Non-dimensional absolute bed temperature. 
Miluid Viscosity . 

Biuwa Kinematic viscosity at minimum fluidization. 
SoOnstant - 

Fluid density. 

Density of packet. 


miata Gensity at minimum fluidization. 


ilies: 





Particle density. 
Ne tiemu cdawed bed densi ty . 


OMe tmemtey son SsOlId particle. 
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If one considers a fluid flowing through a bed of solid 
Barticles, several forms of today's fluidized bed technology 
are created. When the particles are stationary, the process 
is called flow through a porous medium or fixed bed fluidiza- 
tion. If the fluid velocity is such that the particles 
become entrained by the fluid and are removed from the bed, 
the process is then called pneumatic conveying or entrained 
flow fluidization. At some point between these two conditions, 
the particles in mixture with the fluid are transformed into 
a fluid-like state causing continuous interaction between the 
two mediums. The bed particles in this state do not, however, 
leave the bed container. This state is called fluidized bed 
meuaid2Zation. 

The technique of fluidization (all forms) is used for 
solid-fluid interaction processes whenever high rates of heat 
and mass transfer between the two constituents is required or 
high rates of particle transport to and from the bed is 
necessary. The advantages in the use of fluidization can be 
M@eted as follows: [6,12]! 

(a) The formation of an extremely large area of contact 


between the particles and fluid. This very large surface 





‘Numbers in brackets indicate references. 
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area permits the achievement of high overall rates of heat 
and mass transfer between the particles and the fluid. 

(b) The comparative ease with which the fluidized 
particles can be nandled. 

(c) The reduction of temperature gradients to negligible 
Meepertions throughout the bulk of the bed as a result of the 
high degree of particle mixing that can occur in gas-fluidized 
systems, i.e. a very high effective internal thermal conduc- 
memyLty . 

(d) The occurrence of high rates of heat transfer between 
the fluidized particles and an immersed surface. 

(e) The value of the high thermal inertia of the solids 
in heat transfer Operations witnin the bed as a whole. 

(f) The face that the gas-fluidized systems approximate 
a fluid which has a low vapor pressure even at very nigh tem- 
peratures - aS a thermodynamic fluid it lies roughly between 
gases and liguids. 

Inherent witn the fluidization techniques are also some 
mmPeations: [6,12] 

(a) The process demands the expenditure of power for 
mearaization. 

(b) Some solids cannot be fluidized because of either 
an excessive tendency for particle attrition or agglomera- 
meen tO Occur. The first will produce too many fines by 
degradation of the particles which will be lost from the bed. 


Tne latter will give "agglomerated" particles too large to 


1g) 





fluidize. Also, there are those materials of such irregular 
shape that they can only be handled by spouting bed techniques. 

(c) Operating rates are limited to within the range over 
which the bed may be fluidized. If the fluid flow is too low, 
the bed settles; if, on the other hand, the rate is too high, 
there is an excessive loss of particles from the bed by carry- 
over. 

(d) There is usually a limit to the size of the particles 
which can be fluidized. This limit is approximately 6.5 mm 
fees in) down to a few microns (0.001 in). It is particularly 
dependent onthe fluid viscosity and particle/fluid density 
ratio as well as the fluid flow rate. 

(e) True counter flow operation is unobtainable because 
of the high degree of mixing that usually occurs within a 
fluidized bed. 

(f) The general absence of thermal gradients throughout 
a well mixed fluidized bed is a disadvantage in those opera- 
tions where thermal gradients are required. 

(g) The dynamics of fluidized systems are not sufficiently 
well understood to permit confident scale-up to large scale 
units from small scale test data when reaction conditions are 
Seitical. 

Up until this time the term fluid has been used in describ- 
ing the transport constituent in fluidization. To be more 
Specific there are two fluids which have received widespread 


use in fluidized bed technology, air (gas) and water (liquid). 
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The basic differences are in the dynamics and the mass and 
heat transfer of the solid particles suspended in the gas- 
fluidized or liquid-fluidized bed. This is primarily a result 
of the difference in the particle to fluid density ratios. In 
the gas~fluidized bed the ratio will usually be on the order of 
a tnousand while ina liquid-fluidized bed the particles are 
only several times heavier than the fluid. 

Additionally, the fluid flow characteristics will affect 
the overall mass and heat transfer bed characteristics. In 
the liquid-fluidized system, the fluidization velocities are 
usually of such a magnitude that the fluid inertial terms are 
considered insignificant. This is not the case with a gas- 
fluidized system, especially when large particle sizes are 
used. As the velocities increase the inertial fluid losses 
and form drag become important in determining the bed charac- 
teristics. 

With the smooth manner in which liquid-fluidized beds 
Operate, 1.@. no gross instabilities in flow, the bed is 
sometimes referred to as a homogeneously fluidized bed or 
a particulately fluidized bed. The gas-fluidized beds on 
the other hand usually nave unstable flow above some minimum 
meee Zzation condition, characterized by violent agitation 
Of the particles accompanied by bubbling and channeling of 
the gas. As a result, the gas-fluidized bed is sometimes 


called an aggregate or bubbling fluidized bed. 
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The bulk of recent research has been in the area of 
See=tluidized beds vice liguid-fluidized beds. Whether this 
trend will shift remains to be seen. Nevertheless, this 
investigation will be in the gas-fluidized field and all 
further discussions will be based on available gas-solid 
research data. 

In a gas-fluidized bed the particles are supported above 
agas distributor plate which permits the gas to enter and 
flow upward through the bed. As the gas flow increases, 
there will be a pressure drop across the bed which will 
eventually reach a point when it just balances the downward 
gravitational force of the particle bed. When this point is 
reached, the bed will expand and the particles separate. 
This is referred to as the minimum fluidization point. fFur- 
ther gas flow rate increases will result in bubble formation 
which will rise through the bed and burst at the surface not 
unlike a vapor bubble in a boiling liquid. 

The separation of the particles creates a large particle 
Surface in direct contact with the gas and the bubble flow 
induces vigorous and thorough mixing of the particles. This 
high gas-particle contact area results in the high heat 
transfer rates which have been experienced. Because the 
Peetacie flow takes place in the fluidizing gas, a secondary 
fluid-mixing phenomenon will be induced through the particle 
flow pattern. The extent of this secondary mixing and its 


effect on heat transfer will depend on whether the particles 
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are dispersed in the bed more or less as individual particles 
or clumps of particles. As an individual particle, disper- 
sion by the fluid mixing would be relatively small. As an 
agglomerate of particles, however, the result is a dragging 
of gas pockets through the bed resulting in perturbations to 
the particle flow patterns. Thus, the extent of fluid mixing 
is related to the mode and quality of fluidization. This is 
an area that has been scarcely investigated. Two properties, 
though, that have been established as a result of fluid 
mixing investigations are that: 

(a) Owing to fluid back-mixing, temperature driving 
forces are reduced. 

(b) Analogously to heat transfer, driving forces are 
also reduced with mass transfer. 

Based on previous research the normal particle flow 
pattern in small diameter cylindrical gas-fluidized beds is 
as shown in Figure l. For gas flow rates which are signif- 
icantly greater than the minimum fluidization flow rate the 
solids tend to move downward at the bed walls and the passage 
of the individual particle is readily observed. The particles 
will move down the wall and then suomerge into the body of 
the bed. Around this entry point otner particles will emerge 
from tne bed and travel down the wall until they also merge 
back into the bed. This process continues until the bottom 
of the bed is approached. At some point near the bottom of 


the bed an inactive or dead region will form in the particles. 
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The size of this region will depend on the height of the bed, 
fluid flow rate and bed container geometry. In deep beds 

the flow regions are well defined while in very shallow beds 
the downward flow pattern is much less clearly defined with 
the bed appearing to be more turbulent and in a homogenous 
state of agitation. 

In the literature reviewed, when the term shallow bed was 
used it generally was referring to a particle bed whose depth 
was On the order of one and a half feet or less depending on 
the fluid inlet cross-sectional area. In general, the major 
features distinguishing shallow beds from deep beds can be 
described as: [55] 

(a) Violent coalescence of bubbles taking place tnrough 
the majority of the shallow bed depth. 

(b) The jetting region next to the distributor plate 
occupying a Significant proportion of the shallow bed depth. 

mm) The thermal “entry region” occupying a Significant 
proportion of the shallow bed depth. 

These features would suggest that the particle motion 
Within the shallow bed may well be significantly greater 
than in a deep bed and that the bubble growth will largely 
be completed within a very small proportion of the bed depth. 
Thus the residence times of the particles at the heat trans- 
fer surface would be considerably shorter than in a deep bed. 

In the flow description of deep beds, the bed particles 


will move upward in the center of the bed either in a spouting 
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type pattern or in a multi-channel pattern. It should be 
noted that if the bed vessel is Operated in a configuration 
such that its containment walls are not parallel to the 
fluid distributor plate inlet gas flow direction, then the 
particle flow pattern will be adversely affected. Only a 
few degrees of deviation from this parallel configuration 
will lead to particle flow pattern short-circuiting. If the 
Operating gas velocities are further increased while this 
configuration is present, then slugging will initialize. 

As was previously stated, fluid mixing in gas-fluidized 
Meas will occur as a result of the particle mixing at the 
onset of fluidization. Unfortunately, the gas does not 
circulate in the smooth manner shown in Figure 2. Instead, 
there will be many eddies within the central core of the 
particle flow pattern as well as horizontal components 
radiating out from the core. Additionally, part of the gas 
may rise rapidly in the form of bubbles, rising ahead of the 
slower moving interstitial gas components. Therefore, 
Figure 2 should only be viewed as an ideal flow pattern at 
minimum fluidization for the bed. 

To date, particle flow patterns with the resulting par- 
ticle and fluid mixing criteria have been explored only in 
equipment of relatively small diameter (less than one foot). 
The experimental apparatuses that were used were circular 
in nature with either no external heat source or a submerged 


cylindrical (tube) heat source. The data obtained was 
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Mematly a byproduct of the investigation to determine heat 
transfer coefficients for fluidized beds with various heat 
sources and bed geometries. 

This thesis investigation was conducted to determine 
the particle flow patterns developed in a gas-fluidized bed 
of variable rectangular geometry. Secondly, to determine 
what effect an external plate heat source has on the particle 


and fluid flow patterns developed within the rectangular bed. 
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II. OBJECTIVES AND METHOD OF APPROACH 


A. GENERAL STATEMENT OF THE PROBLEM 

The specific problem may be defined as an investigation 
into the formation of particle flow cells at minimum gas 
fluidization flow rates in a rectangular gas-fluidized bed. 
The primary objective is to determine if multiple cell for- 
mation occurs as the bed width to height (W/L) rao 
increases. The secondary objective is to determine if the 
presence of plate heaters at the end walls of the particle 
bed will alter the particle flow cell formation patterns 
when energized and to what extent they affect the fluid 
Mixing (heat transfer) flow patterns within the particle 


LOW . 


Bee METHOD OF APPROACH 

The methods used for this investigation were experimental 
and analytical. Using available research literature a deter- 
Mination was made of the minimum fluidizing velocity (u¢) 
For the particular fluidized bed configuration. This value 
Of Ung Was based on recorded and analytical pressure drops 
across the distributor plate and bed in addition to the 
expanded bed heights. Visual confirmation of the minimum 


fluidizing conditions was then compared with the analytical 


estimates of the required pressure drops, bed height and air 
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velocities. Once an initial particle flow cell was estab- 
lished, the bed geometry and/or activation of the plate 
heaters would be used in an attempt to alter its configura- 
tion. A grid of variable depth copper-constantan thermo- 
couples was used to map the thermal distribution within the 
particle flow cell to determine the fluid mixing properties 
moeene cells. 

This thesis will not be involved with the heat transfer 
coefficient determination so prevalent in recent experimenta- 
tion; rather it will focus only on how the particle flow cells 
relate to the bed geometry chosen and how the resulting fluid 
mixing affects the heat transfer within the bed. It should 
be noted that direct measurement of the solid particle temper- 
atures 1S very unreliable and that the temperatures recorded 
fee be considered to be that of the air at a particular 
point within the grid of the bed. Because of the difference 
in the particle to fluid density, however, the two will be 
considered to be at equilibrium at the point measured for 


this experiment. 
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TTI. SUMMARY OF PREVIOUS RESEARCH IN 
Gla TUS DIZED BED TECHNOLOGY 
A. GENERAL DESCRIPTION OF GAS-FLUIDIZED BED BEHAVIOR 
1. The Phenomenon of Fluidization 

When the bed configuration is initially filled with 
particles, an initial voidage (€ 4) will result between the 
particles. As the fluidizing fluid, a gas in this case, is 
initiated at some superficial velocity (u 4) there will be no 
movement of the bed. However, a pressure drop (Ap 3) will 
occur and gradually increase at a linear rate. The magnitude 
of this pressure drop will be determined by the gas flow rate 
and the characteristics of the bed. 

Eventually a condition will be reached when the gas 
flow rate (a) has increased to a point at which expansion 
of the bed initiates. The voidage will increase to some Eo 


and the pressure drop can be calculated from the equation 


ne eee) GO Sule) ) (alg 


(Note - The units to be used with all text equations can be 

found in Appendix D.) While the bed is expanding, there 

Still will have not been any particle motion taking place. 
As the gas flow rate is increased further, a point 


will next be reached when the voidage becomes ¢€ and the 


mf 


particles are at the point of fluidization onset. The 
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velocity is called the minimum fluidization velocity (ug) 
p@eeis Critical to most bed behavior calculations. Further 


gas velocity increases (u,) will fluidize the bed causing 


Q 
further voidage increases (Eo) until a terminal velocity 
(u,) is reached as defined by the bed container geometry or 
particle entrainment takes place. By the time uy has been 
reached the particle motion will be intense although the 
pressure drop across the bed will remain at the level of 
minimum fluidization inception. This condition is referred 
to as dense-pnase turbulent fluidization at its inception 
(U5) but degenerates to a dilute-phase when the gas rate 1s 
Seen that the particles become entrained in the gas. On the 
average the ratio of the terminal velocity (u,) to the 
ieeamum fluidization velocity (ue) is about seventy (70). 
The terminal velocity is also a function of the particle 
diameter, shape and density as well as the properties of the 
medi zing fluid. 

Parent et al. [52] has characterized three types of 
behavior with the gas-solid mixtures in the dense-phase 
MmmeadiLzation region. Under certain conditions, usually for 
small particle sizes, the gas tends to channel through the 
bed without cauSing movement of the bed. Under other condi- 
tions, such as large bed heights, large particle sizes and 
high gas flow rates, slugging takes place which is analogous 


to the turbulence observed ina boiling liguid. The third 


type of behavior is described as the mechanically smooth 
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Fluidization. In this case the particles move along twisting, 
ever-changing paths. Occasionallly gas bubbles appear in the 
bed but there 1s no rough bumping. It is this case that the 
pressure drop across tne ded 1s nearly equal to the suspension 
head opposing the flow (drag). A word of caution is warranted, 
however, in that in research conducted by Lewis et al. [42] 
the conclusion that the pressure drop is essentially equal to 
the weight of the particle bed would appear to be satisfactory 
only for units with small L/D ratios. It appears to be in 
considerable error for large L/D values or for beds having 
large static effects present. 

Here again several researchers differ in their opin- 
ion as to the cause of this phenomenon. While Lewis et al. 
[42] explained it as an indication of frictional drag of the 
solid particles on the walls of the unit, Toomey and Johnstone 
[62] attribute it to kinetic energy losses caused by the colli- 
sion between the particles and Wilhelm and Kwauk [76] con- 
cluded it was caused by forces which are electrostatic and/or 
Fluid dynamic in origin. Since this thesis is operating in 
the shallow bed region it 1S not anticipated that this excess 
pressure drop will be of concern. 

2. Voidage at Minimum Fluidization 

The first bed characteristics that must be determined 
for an understanding of gas-~-fluidized bed operation is that 
of the bed voidage at minimum fluidization (e ¢)- By defini- 


eron E mt is that particle separation to which the bed must 
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expand before particle motion, as induced by fluid flow, can 
initialize. There are two methods whicn can be used to deter- 
feme this factor. 

The first method is based on a collection of data from 
several researchers that has been graphically combined by Leva 
Mioweand reproduced here in Figure 3. For a given particle 
diameter (D,,) and composition a value for the minimum fluidiza- 
tion voidage (€ £! can be obtained. The obvious drawback for 
this method is that it depends solely on the material chosen 
for the particles being one that has had previous experimenta- 
tion accomplished. 

The second method is used with those beds for which 
particle te data is not available. It requires that the bed 
be subjected to the fluidizing gas until particle motion 
initializes, at which time the bed height (Le) 1s recorded. 
The minimum fluidization voidage (€ Ve) can then be calculated 


from the equation 


E = | - —-————  ———_ —__ (eoneeee) 
mf 
The accuracy of this estimate will depend on the mass charac- 
teristics of the bed particles, i.e. sphericity, whether solid 
Seenellow, degree of particle porosity, etc. Particles of 
high porosity are usually unsuitable for fluidized bed use 


as they have a tendency to break up with increased fluid flow. 
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Figure 3. Values of €pr vs. particle diameter. 
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Once a value for the minimum fluidization voidage 
has been estimated, additional bed behavior characteristics 
can be predicted and the bed's suitability for the intended 
process determined. 

meee Gessure Drop 

Another excellent indication of fluidization behavior 
of the particle material may be obtained from the pressure 
drop-flow relationship of the expanded bed [36,40,64]. Ideally, 
the pressure drop across the bed should increase linearly until 
the voidage has reached a maximum value Ce.) at which time the 
pressure differential would level off for the remaining phases 
of fluidization. See Figure 4. 

Unfortunately, the idealization never occurs as a 
result of three major factors; namely channeling, slugging 
and distributor design. With channeling a pressure drop 
occurs (Figure 5) when 2a is reached. The magnitude of this 
drop will depend on the extent of the channeling phenomenon 
but eventually a minimum drop point is reached and the pres- 
Sure differential then increases to some new lower steady 
State value. 

The channeling can be either "through channeling" in 
which it extends completely through the bed or "intermittent 
channeling" where it extends tnrough part of the bed. For- 
tunately, as the gas flow rate is increased the tendency for 
the bed to have channeling decreases. Some factors [38] 


Which tend to promote channeling are poor distributor design, 
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Figure 4. Ideal pressure-drop-flow diagram. 
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Figure 5S. Pressure-drop-flow diagram for channeling. 
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the shape of the bed particles (less sphericity the greater 
the channeling), moisture in the bed and the rate of fluid 
flow through beds of small particles. If the bed is moist, 
the particles will tend to clump together, and the gas flow 
will channel between the clumps. At higher flow rates the 
increased agitation of the bed apparently destroys the 
channels as soon as they are formed. 

With slugging a condition exists in which the bubbles 
of gas combine to a size approaching the order of magnitude 
of the diameter of the confining vessel. As a result, the 
particle layers or slugs of granular particles will move 
upward through the bed like a piston until they disintegrate 
at the bed top. When this happens, the pressure drop will 
increase above the predicted steady state value (Ap), see 
Figure 6. If the flow rate is increased further, the pres- 
Sure drop rise continues and the bed will slug even more 
severely. 

Bie teninad factor, namely distributor design, was 
well investigated by Grohse [27] who stated that for the mul- 
meorirtice plate the point of incipient fluidization is not 
well defined. The gradual pressure-drop increase up to the 
theoretically demanded suggests that the bottom portion of 
the bed is not fluidized to the same extent as the upper 
region. With the screen distributor, a certain amount of 
Channeling probably occurred. There is the possibility that 


the holes in the screen were locally blocked; this could 
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create a condition favorable for channeling. With the porous 
plate the point of incipient fluidization is well defined 
and there is no tendency for channeling. See Figure 7. 

Pacem cne point Of minimum fluidization is the key 
Starting point on fluidization studies, it seems extremely 
important to be able to predict that point both theoretically 
and through experimental results. Therefore, the distributor 
design's effect on the pressure drop at Uf becomes a very 
important factor. According to Kunii and Levenspiel [34], 
distributors should have a sufficient pressure drop to 
achieve equal flow through all openings. As a result, the 
pressure drop across the distributor plate by direct measure- 
ment should be at least 10% of the pressure drop of the bed 
to ensure effective fluidization. Theoretically the pressure 
drop at the minimum fluidization velocity (AD «¢) can be 


found by substituting Lf and ¢€ for the values of L, and 


mf 
an mamequation (3.1). 
fee Minimum Fluidization Velocity 
Mie sdollpey tO predict the point Of initial fluidi- 
Maeton reliably is of basic importance in virtually all 
fluidized-process studies and designs. By operating at mass 
velocities which are slightly higher than the minimum fluid- 


ization mass velocity (G auata bionwand Clutrirakien 


we 
(stratification) losses will be minimized. The by-product 
of the low flow rates is then a low particle velocity in the 


bed and corresvondingly a relatively low bed wall heat transfer 


Boer ricient (hi). 
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The general equation [36] for the minimum fluidization 


mass flow rate 1s 


r 1.82 ~ 
688 D., { Oe ( 0. O¢ 
G a mm mm mm mn ’ Cees. 
mf 1 098 8 ( 


with Oe the average particle size, 


D = - (34) 





Equation (3.3) can only be used, however, when the Reynolds 


number 


0.03445 D., Gat 


oe (S95) 
u 


is less than ten. Should the value of Re wf be greater than 
ten, then the alternative is to use the nomograph [72] of 

Figure 8, with the calculated value of Dor u and p¢lo, - Ae) 
to recalculate the value of Ge’ If, however, the value of 
mo, u or Pelo, = Pg) does not appear on the nomograph, then 


P 


go to Figure 9 [72] and with the calculated value of Re wf! 


from equation (3.5), determine the correction factor to be 
applied to the value of Ga Ga culated from equation (3.3). 


The new value of Gt 1s then calculated from the equation 


( oe ) pee ae Gut Pe omeaelon ( Gat 


) (3:46)) 


calculated’ 
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It should be mentioned that there are numerous corre- 
lations proposed for the minimum fluidization mass velocity 
(G_¢) that are applicable to the experimental results of 
various researchers. These predictions are not in very close 
agreement and are a function of different apparatus construc- 
tions, fluids used, particles used, etc. 

One variation considered worthy of note is a result 
of the work accomplished by Wen and Yu [73]. They demonstrated 
that when particles of more than two different sizes are fluid- 
ized, they either separate completely, forming distinct layers 
Seeintermix forming one layer of rather uniformly distributed 
particles. For a bed containing two different particle sizes, 
if the particle diameter ratio is less than 1.3, the particles 
were observed to intermix. If, however, the particle size 
ratio was greater than 1.3, the smaller particles fluidized 
separately above the larger particles. When this occurs, 
the voidage must be obtained for each layer and special pro- 
cedures, delineated in Reference [73], used to compute the 
fluid mass velocity. 

5. Bed Expansion 

This area of bed behavior is important in that it 
resolves the minimum bed chamber height that must be designed 
mBomcenable fluidization to occur. Unfortunately, although 
gas-fluidization processes are more prevalent than liquid- 
fluidization processes, the mechanics of the gas-fluidization 


process is not as well understood in comparison to tne liquid 
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process. It has been found that the liquid systems will ex- 
pand smoothly almost to a voidage of 1.0. With a gas system, 
however, somewhere along the course of the expansion, a point 
of mechanical instability is reached which degrades the expan- 
Sion process. The location of this trip point or the reason- 
ing as to why it occurs has so far eluded researchers. There 
is some speculation, however, that it may be caused by the 
formation of aggregates of particles in the bed and that it 
1s somehow related to the fluid velocity. 

Based on analytic and experimental data, Leva et al. 
[36,40] formulated a simple relationship between the fluid 
mass velocity and bed particle voidage that enabled the 
generation of a prediction model for bed expansion. Some 
examples of this gas-solid expansion data can be found in 
Figure 10. What is shown is that at higher flow regions the 
data will align along curves rather than linear lines and at 
intermediate flow regions the reverse is true. Note that as 
tne vOidage increases, (1 - ¢€)/e* will decrease. Also, at 
elevated flow rates, the slopes (m) become less negative in- 
dicating that aggregation becomes more prevalent while the 
linear lines indicate little or no aggregation (intermediate 
BbOw region). It is also noted that small-particle beds may 
be expanded over substantially wider flow ranges than larger- 
particle beds before aggregation conditions will set in. 
This 1S indicated by the fact that, with smaller particles, 


departure from the straight line occurs at higher values of 


the reduced mass velocity than with larger particles. 
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Figure 10. Gas-particle expansion data. 
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As far as practical methods of correlation are con- 
cerned, the data is most cases may be represented by straight 
lines. This will produce the relationship between slopes (m) 
and particle size (D,,) indicated in FPigure ll [36,40]. This 
correlation is based on virtually all of the available bed- 
expansion literature data to date. 

To calculate a predicted expanded-bed height of the 
gas-fluidized particles, two basic methods may be used. The 
choice of method will depend on the bed-particle diameter. 
However, a Maximum mean particle diameter (D.,) Orrapout 0-015 
in is recommended, otherwise the fluidization efficiency will 
be so low as to promote slugging. Fluidization efficiency 


(n) is defined by the equation 
es ’ CBee?) 


where G. and en are respectively defined as the fluidizing 
gas mass velocity flowing through a bed at a given instant in 
time and the gas mass velocity required to initially expand 
the bed. Once slugging sets in the bed expansion height 
fate begin to fluctuate as the gas bubbles expand and collapse. 
When this observed condition occurs no correlation to date 
has been found to give a predictable bed height. 
a. Method One 
(1) With D, known, m may be found from Figure 1l. 


(2) Determine E nf PrOonereaure 3° Or Equation (3.2). 
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(3) Calculate Gat uSing equations (3.3) and/or 
meeo) Or Figure 8. 

(4) Draw fluidization line on log-log paper with 
slope m, see Figure 12 [40]. 


(5) For the cnosen value of G find the corre- 


‘ee 
Sponding value of Eee 
(6) Solve the voidage equation, equation (3.2), 
Or Le Dy SUDSEIEUtInNg L and E - BONG Lit and emf’ 


b. Method Two 

(1) Calculate Ge uSing equations (3.3) and/or 
mo) Or Figure 8. 

(2) Calculate the reduced mass velocity, G,/G_-- 

(3) Go to Figure 13 [36] with the reduced mass 
velocity and Bs to find the fluidization efficiency (n). 

(4) Go to Figure 14 [36] with n and D., to evaluate 
the bed expansion ratio, R. 

(5) Calculate the static bed height Le from 
Semeation (3.2). 


ey) uibulss slowly, Ls by R to get the expanded bed 


1 

height, Lee 
Once the expanded bed height has been calculated 

there still remains one check to be performed on its validity. 

This is a check on the extent of height fluctuations that may 

Occur for the flow rate chosen. It has been found that under 


conditions of high gas flow rates that an oscillation of the 


expanded bed height can occur. This concept has been termed 
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Figure 13. Fluitdization efficiency vs. reduced mass velocity. 
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the fluctuation ratio, r, the ratio of the highest to lowest 
levels which the bed top occupies for a particular gas flow 


rate. The value of r can be determined from the equation 


mo Se > (oe 
r=>exp --——_—_—_——— (3.8) 
G 
mf 


where slope m*~ is related to the particle diameter as shown 
in Figure 15 and is based on data of Lewis et al. [42]. The 
Seoser the fluctuation ratio to 1.0 the more accurate the 


expanded bed height estimate. 


B. HEAT TRANSFER PHENOMENON MODELS 

There is little agreement in the reported experimental 
values of the heat transfer coefficients in the literature. 
The reasons can best be summarized as follows: 1) a variety 
of contacting patterns between particles and the fluid re- 
sulted in different models from which to base the average 
heat transfer coefficient; 2) the temperature measurements 
taken in the beds by thermocouples were interpreted differ- 
ently (usually the readings were related to the gas temper- 
ature in some manner); or 3) that direct measurement of 
sOlid particle temperatures are unreliable with present 
state-of-tne-art methods, especially for the smaller size 
particles (less than 100 microns). 

As a result, due to the generally low accuracy of in-bed 


temperature measurements and the oversimplifications in the 
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Figure 15. Slope m’ in relation to particle diameter. 
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flow models, heat transfer coefficients so far reported have 
varied widely. Although there is no complete theory that 
can predict the heat transfer data for fluidized beds, there 
are several theories that correlate the published data. Each 
resulting heat transfer model has some experimental or design 
Situation for which it 1S appropriate. What follows is an 
Overview Of these various models all of which are based on 
cylindrical bed geometries. 

1. Steady State Conduction Across the Gas Film 

In this model heat is conducted through the gas 
boundary layer near the heat exchange surface. In order to 
explain the high values of the heat transfer rates obtained 
experimentally, the gas layer is assumed to be scoured by 
sOlids moving along the heat exchange surface, thus decreas- 
ing the boundary layer thickness near the wall. This model 
was developed by Leva et al. [37], Dow and Jakob [13], and 
Levenspiel and Walton [41]. 

Dow and Jakob [13] studied the gas-fluidization mix- 
tures and classified them as either a dense or lean phase 
mixture based on the gas velocity. The dense phase occurred 
at gas velocities close to Une while the lean phase was close 
to u,. Their model applied to the mechanically smooth dense 
phase fluidization although they did observe bubble formation 
and slugging in the dense phase region of fluidization. 

The model's motion, Figure 16, is primarily upward 


at the center of the bed and downward along the wall. 
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Temperature measurements snowed the entire bed was at a con- 
stant temperature except for a small region at the bottom 

and a narrow layer near the wall. Heat is tranferred radially 
tnrough a thin air film to the particles and air moving down- 
ward outside this film. The particles carry the heat to the 
bottom of the bed where thermal equilibrium is attained almost 
instantaneously with the incoming cold air. The bombardment 
of small particles prevents the formation of a laminar bound- 
ary layer, and leaves a very thin laminar sublayer and a 
thicker turbulent outer layer. 

Tne heat transfer rate from the surface to the fluid- 
ized bed was found to be independent of the thermal conduc- 
tivity of the solid particles, and tne large isothermal 
region of the fluidized bed could be destroyed if the down- 
ward flow of the small particles was short-circuited into 
tne main flow before they reached the bottom, i.e. under 
slugging conditions. 

Levenspiel and Walton [41] presented a similar model 
where the resistance to the heat flow is due to the laminar 
Gas layer which is destroyed by solid particles passing 
through it. Thus the average thickness of the laminar layer 
is much less than in an empty tube. They assumed that the 
particles are of uniform diameter ea Stationary, and 
arranged in equally spaced horizontal layers. The boundary 


layer formation at the wall 1s similar to that of a flat 
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plate, see Figure 17. At points of contact with the station- 
ary particles, the boundary layer is destroyed and starts 
once again. 

Unfortunately, the two models do not take into account 
the influence of the solid particles on the heat transferred, 
so they must be considered of limited application. 


2. Unsteady Heat Conduction by Single Particles in 
Direct Contact with the Heat Exchange Surface 


With the advent of photographic studies of the bound- 
ary wall during fluidization, it was shown that the particles 
tended to associate in groups which moved slowly along the 
wall or remained in contact with the wall a certain length of 
time. Those observations led Botterill et al. [6] to propose 
that heat was transferred from the vertical surface to the 
adjacent layer of particles arriving together at the wall. 

If the residence time of the particle near the wall was long, 
heat would penetrate further into the bulk of the particles 
and if the residence time was short, heat exchange calcula- 
tions would be limited to only one particle, the nearest to 
the surface. They confined their theoretical considerations 
to an isolated particle lying in contact with the wall. 

After much experimentation this model was extended 
to cover different particle sizes and materials and a wide 
range of residence times for the particles at the surface. 
The results showed the validity of the model for short resi- 


dence times with the observed effect that the neat transfer 
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Figure 16. Heat transfer mechanism and particle motion. 
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coefficient 1s a decreasing function of tne residence time. 
meen, this model can only be applied to mechanically smooth 
dense phase fluidization as the model's validity breaks down 
in the lean phase or when bubbling/slugging occurs. 

3. Unsteady Transfer of Heat to "Packets" of Particles 

Which are Renewed by Violent Disturbances in the 

Core of the Fluidized Bed 

In an attempt to present a mechanism for the heat 
transfer in bubbling, dense phase beds, Mickley and Fairbanks 
[47] proposed a model of a small group of particles moving as 
individual packets. They based their modelon the fact that 
tne voidage of the dense phase bubbling beds was found to 
be close to the voidage of quiescent beds. 

In the dense phase fluidized bed, each particle was 
expected to be in contact with several neighbors most of the 
time. The packets are not considered permanent; they Nave a 
finite life expectancy; and their voidage, density, heat con- 
ductivity, and heat capacity are assumed to be the same as 
those of the quiescent bed. See Figure 18. 

This model assumes that as a packet at bed tempera- 
ture Ty comes into contact with a flat surface at a tempera- 
ture TO that unsteady heat transfer will start upon contact. 
Meeeher, if An is the contact area between the packet and the 
Wall, and the packet is homogenous, then the heat transferred 


after some time t will be given by the equation 


ee) 





miaere the local instantaneous heat transfer coefficient is 


given by the equation 


lee = — e eSpethiG)), 


Mickley and Fairbanks [47] also presented two devia- 
tions from their original model which covered the phenomenons 
of slug flow and side mixing. The first, on slug flow, was 
for fluidized beds with very low gas flow rates and no turbu- 
lent mixing. In this model, Figure 19, the particles were 
to move downward along the outside wall of the bed. The flow 
is not smooth but is spasmodic owing to the intermittent 
passage of ascending gas bubbles near by in the core of the 
bed. However, assuming that the particles moved downward at 
a constant velocity (ul), the age of all the packets at a 
distance & from the bed top was t = &/, and the resulting 
average heat transfer coefficient for a heater of length L 


H 


could be found from the equation 


= : Coe 1) 


The second model deviation, on side mixing, was for 
beds with large heat transfer surfaces that were highly 
turbulent. To represent this case, Figure 20, they assumed 
that as the particle packet is moving downward some of it 


is exchanged with some other packet at temperature Th? brought 
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in sideways from the core of the bed. The average number of 
replacement packets at the wall per unit time resulting from 
Side mixing was then defined as s, resulting in the local heat 


transfer coefficient equation 


2 0.5 tow 2) 
DL ( Kin om mn Ss) 


Tt should be noted that for these penetration type 
models the residence time of the particle is of fundamental 
importance and has generated a considerable amount of research 
both in its determination and its effect on the various models 
Sie this section. 

Patel [53], in a further spin-off of Mickley and 
Fairbank's Original models presented two surface renewal 
models that can be characterized as follows. 

In model one, Figure 21, a packet initially at the 
bed temperature TS arrives at the wall (temperature Te) at 
time t = 0. The packet is assumed to have the same proper- 
ties as the bed at minimum fluidization. It receives heat 
from the wall through a contact resistance and returns to 
the bed after some length of time. During that residence 
time, heat penetrates a distance x; beyond this distance 
the temperature is constant and equal to the bed tempera- 
ture Th: 

In model two, Figure 22, the fluidized particles 
are assumed to be spherical and of the same diameter (d} 


A particle at the bed temperature T, arrives at the heating 


b 


6 3 








Tiree Mais AON 


mroere <4. Downtiow of particles with side mixing. 


Packet 
Wall (Ty) apes Gt.) 
Canwac t mit ial 
Resistance Temperature 
Cie i 





micume ci ratelevodel. I. 


64 





surface (temperature qT.) and receives heat from the fluid 
Mayacent to the wall. This fluid is assumed to be stagnant 
and at wall temperature De While the particle is at the 
Surrace it also loses heat by conduction to a packet of par- 
ticles of thickness x situated between the wall particle 
and the bulk of the bed. The packet is assumed to have the 
Same properties as the bed at minimum fluidization conditions. 
Patel by no means had the only variation on Mickley 
and Fairbank's packet theory. Numerous other variations 
exist with different Dee values being correlated by each 
researcher. Appendices A through C give some of the result- 
mig correlations. 
meeoceaay Convection Through the Emulsion Layer 
VeaneHeerdeneect at. [OOP Observed that the heat capa- 
city of the particles per unit volume is of an order of mag- 
nitude of thousands times greater than that of the gas and 
that the mean particle velocity is much lower than the gas 
velocity. Therefore, the largest portion of the heat will be 
transferred by the moving particles. They assumed that the 
gas temperature was that of the particles and as a result 
that the gas only provided a suspension for the particles. 
They observed negligible radial motion of the bed 
particles and therefore concluded that the radial heat trans- 
fer was a function of the tnermal conductivity of the sus- 
pension. Experimentally, they verified that the heat transfer 
coefficient is large for short heat exchange surfaces and vice 


versa. 
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Wicke and Fetting [75] furthered van Heerden's work 
and proposed a model, Figure 23, for those observations. [In 
their model, heat (q_) from an exchange surface was first 
transferred by conduction through a gas layer whose thickness 
was SG: This heat was then divided into two components: 
namely heat (q,) taken by the particles flowing parallel to 
the surface in a second zone of emulsion of thickness Sor and 
heat (q,) transferred into the core portion of the bed by 
interchange of the particles. In their calculations, how- 
ever, they neglected the term qe: 

The last model to be presented is actually a combina- 
tion of several heat transfer mechanisms proposed in earlier 
models. 

5. Kunii and Levenspiel Model [34] 

In this model, Figure 24, four mechanisms operate 
together. 

a. A film of gas coats the surface. Its thickness 


(6.) is large or small depending on whether a bubble of air 


G 
is near the surface or the emulsion is uniform and close to 
the surface. 

b. Some solid particles are in direct contact with 
tne heat exchange surface. 

c. There is a layer of emulsion with thickness 5 
Which flows along the wall. 
d. Part of the emulsion layer is replaced occasion- 


ally by fresh emulsion coming from the core of the bed or by 


bupbles rising along the wall. 
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Figure 23. Wicke and Fetting Model. 
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Depending upon the fluidization conditions and the 
position of the heat exchange surface in the bed, emphasis 
can be placed upon any combination of the above four 


mechanisms. 


C. FLUIDIZED BED VARIABLES THAT AFFECT HEAT TRANSFER 

In addition to the large number of models that have been 
proposed by various researchers to explain their experimental 
results in neat transfer, there are a number of variables 
which have been found to affect the heat transfer rate. In 
general, the factors which will control the heat transfer 
characteristics of the bed will be those variables which 
determine the fluidizing properties of the bed. In this 
section an attempt will be made to list some of the more 
important factors and their effect on heat transfer. An 
excellent indepth review of general fluidized bed character- 
istics was researched by Gutfinger and Abauf [28]. Addition- 
ally, Appendices A through C present a list of nondimensional 
correlations reported by different investigators while trying 
to correlate their experimental findings. 

ie Properties of the Fluid 

a. Density (p ¢) and Specific Heat ‘Soa OG. Ehewr laud 
The dependence of the heat transfer coefficient 

on the specific heat of the gas has not been clearly deter- 
Mined as the existing data shows no common trend. However, 


Gelperin and Einstein [18] have predicted that there would 
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be an increase in the heat transfer rate with an increase in 
Cor Sp the product ge oe when the fluid was at high pressures 
and velocities. 

Bb: Gas Viscosity (u) 

The majority of researchers agree that with an 
increase in the gas viscosity the wall heat transfer coeffi- 
Client will decrease. 

c. Thermal Conductivity of the Gas (kK ¢) 

This is the fluid property that has the most pro- 
nounced effect on the heat transfer coefficient. Depending 
On which researcher is consulted the value of he will increase 


Peoportional to k. to the power of 0.5 - 0.67. 


if 
2. Properties of the Bed Particles 
a. Diameter (D.,) and Shape of the Particle (>.) 

A significant number of researchers [37,47] have 
concluded that the heat transfer coefficient decreases when 
large diameter particles are used. However, the results of 
Sarkit et al. [60] as reported by Gelperin and Einstein [18] 
show that for the laminar flow region the heat transfer 
coefficient varies inversely with particle diameter and for 
the turbulent flow region it has a direct variation. 

The values for the heat transfer coefficient tend 
to increase as the sphericity (>.) increases (rounder and 
smoother particles). 

b. Density of the Particle (9 .) 
Experimentation has shown that the wall heat 


transfer coefficient increases with the increasing density 
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of the solid particles. However, this increase slows in 
its magnitude as the fluid flow becomes more turbulent. 
G. Specific Heat of the Particle (ea! 

The heat transfer coefficient increases with an 
increase in the specific heat of the particle. 

d. Thermal Conductivity of the Particle (k .) 

Based on measured experimental results over a 
wide range of particle thermal conductivities it has been 
determined that the effect of the thermal conductivity of 
the particle on the heat transfer coefficient is slight, so 
much so that it is usually neglected as a factor of concern. 

fw Luidization Conditions 
a. Gas Mass Velocity (G) 

The increase in the heat transfer coefficient is 
BeEOportional to Gc” where n will vary according to the experi- 
ment being conducted. Its range has been shown to be quite 
large. About the only facts that have been agreed upon by 
researchers for this variable is that of the general shape 
of the gas mass velocity versus heat transfer coefficient 
curve [72], Figure 25, and that there is some maximum value 
wor Ds which exists at some gas mass velocity Gopt° Tables 
listing the nondimensional correlations for both the rising 
and falling sides of the curve, for the maximum heat transfer 
coefficients, and for the special gas flow rates when the 
h is attained are presented by Einstein and Gelperin 


wmax 
[14]. 
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Figure 23. Heat transfer coefficient vs. mass velocity. 
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b. Bed Porosity or Voidage (ec) 

Based.on film theory the neat transfer coeffi- 
cient will vary inversely with some function of the bed 
porosity. That is, as the particle aggregate at the wall 
decreases so will the heat transfer rate. However, there 
has been some experimental data which has demonstrated the 
opposite effect. 

4. Geometric Variables 
a. Bed Diameter (D) 

It has not been demonstrated that there exists 
a dependency of the heat transfer coefficient on the bed 
diameter. 

b. Diameter of an Immersed Heater (D,,) 

Researchers are in general agreement that the 
wall heat transfer coefficient increases with a decrease in 
the immersed heater (tube) diameter, and that the heat trans- 
fer rate diminishes as the diameter of the heater is increased. 
Tnese effects are valid up to a particle diameter of 10 mm 
PONS9 in) after which ie is independent of the heater 
diameter. 

c. Fluidized Bed Height (L) 

There has been no dependence shown between the 

heat transfer coefficient and the fluidized bed height. 
ad. Length of Heat Exchange Surface (L,,) 
This is a variable that has had several differ- 


ent results on its effect on the heat transfer coefficient 
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all of which have been substantiated by additional research. 
Basically, one research group has found that the heat trans- 
fer coefficient decreases with an increase in the length of 
the heat exchange surface. Another research group concluded 
Emat if L/D was greater than 7 for an external heat transfer 
surface, then the heat transfer coefficient was independent 
of heater length and that internal heat exchange surfaces had 
no observed effect on the heat transfer coefficient. And yet 
a third research group found that if the heat exchange surface 
was located above the active section of the bed (mixing 
region), then the heat transfer coefficient would be inde- 
pendent of the heater's length. 

e. Heat Excnange Surface Roughness 

It has been determined [31,68] that the heat trans- 
fer coefficient is sensitive not only to the type and orient- 
ation of the roughness but also to the size of the particles 
relative to the size of the roughness. When the particles 
are of a large size relative to the roughness, the gas con- 
duction paths will be longer and the heat transfer will be 
lower aS compared to a smooth tube or surface. See Figure 
26a and 26b. 

When particle and roughness sizes are similar, 
the multiple contact should reduce the mean gas conduction 
paths and gives rise to a considerable increase in the heat 
transfer, Figure 26c. However, possibly due to particle 


Capture and the consequent increase in particle residence 
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times, the increase in the maximum heat transfer coefficient 
for vertical grooves, 1S of the order of 50% for the smaller 
particles and 18% for the larger particles. The particle 
capture and inability of the bubbles to replenish particles 
in the case of horizontal grooves undoubtedly contributes to 
this difference. 

Where the size of the particles is much smaller 
relative to the roughness size, particle capture is not a 
serious problem, but there is some improvement in particle- 
surface contact as well as a larger effective area for heat 
mrensfer, Figure 26d. 

5. Miscellaneous Variables 
a. Entrance Effects 

Entrance effects ina flow situation cannot be 
ignored or avoided. Control or limiting entrance effects can 
only be accomplished by careful design of the gas distributor, 
which in turn, will influence the configuration of the active 
(mixing) region of the fluidized bed. 

b. Temperature and Pressure of the Fluidized Bed 

The heat transfer coefficient is independent of 
the pressure of the fluidized bed or fluidizing fluid. The 
temperature of the fluid, however, is another matter. In 
general it can be assumed that the heat transfer coefficient 
is independent of the temperature of the bed or fluid. 
Studies by Zabrodsky [79] and Yoshida et al. [77], however, 


indicate that the radiation component of the overall heat 
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transfer coefficient must be taken into account when the bed 
Eemperature exceeds 900-1300 K. In turn, the magnitude of 
the radiation component effect will be dependent on the bed 


container geometry. 


D. IMMERSED BODIES 

Experimental work conducted in fluidized beds can be sub- 
divided into two distinct categories. The first deals with 
all experimental work that uses an external heat transfer 
surface, 1.e. the wall of the fluidized bed container itself 
is used as the heat exchange surface. The second category 
uses an internal heat transfer surface. In this case, bodies 
of different shapes (sphere, cylinder, plate, coil, etc.) 
which are immersed in the bulk of the fluidized bed become 
the heat transfer surface. 

While the previously discussed variables and character- 
istic bed behavior are equally germane to either category, 
the use of immersed heating surfaces will have additional 
factors that will need to be taken into account during a 
design process. The location of the heat source with respect 
to the gas distributor plate and its plane of orientation 
Will have a significant effect on the heat transfer coeffi- 
cient. The results of the variations are sometimes contra- 
dictory and depend strongly on the local hydrodynamic 
conditions in each case. 

Limited experimentation has been conducted in the use of 
Spheres, flat plates and small cylinders with the data 
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obtained being sketcny at best. The bulk of the immersed 
body experimentation has been with tubes for which additional 
variations in heat transfer tube design (size, shape, spacing, 
gap, pitch and material) will complicate the general heat 
transfer mechanism. An excellent indepth coverage of this 
topic was conducted by Saxena et al. [61] and is summarized 
as follows. 
1. Horizontal Tubes 

The total heat transfer coefficient cae) for single 
cylindrical tubes increases with factors that increase the 
gas thermal conductivity near the tube wall (i.e. temperature 
at the tube surface) and reduces the particle residence time 
(i.e. decrease in tube diameter and increase in fluidization 
velocity). However, at higher gas velocities the decrease 
in bed density reduces the heat absorption capacity of the 
fluidized particles, resulting in a decrease in the heat 
transfer coefficient (hi). The heat transfer coefficient 
also increases with decreasing particle diameter eo 

The heat transfer coefficient is highest for tubes 
with a lenticular cross section for dense moving packed beds. 
Gas pockets are not created in the flow wake (downstream 
Side) of a lenticular tube, and the particles are not stag- 
nant on the upstream side of the tube. 

The presence of surrounding tubes in an array affects 
the total heat transfer coefficient for an individual tube 


by resisting particle movement and by altering the particle 
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temperature from that of the bed bulk, which reduces the tem- 
perature driving force. This is at least one of the reasons 
why the heat transfer coefficients for the lower tubes are 
S@reater than for the upper tubes, and a close spacing of the 
tubes results in a decrease of the heat transfer coefficient. 

Fins will increase the rate of heat transfer to the 
tubes. However, their effectiveness increases, reaches a 
maximum, and then decreases with decreasing spacing distance 
and increasing fin length [54]. 

An increase in tube-to-tube spacing in a bundle of 
short-finned tubes, which reduces the resistance to particle 
movement between the tubes, increases tne total heat transfer 
coefficient. However, an increase in tube-to-tube spacing 
for a bundle of long-finned tubes does not noticeably in- 
Crease the heat transfer rate, implying that the resistance 
to particle movement between fins predominates. 

2. Slanted Tubes 

There does not appear to be any particular advantage 
in slanting the heat transfer tubes. Minimum heat transfer 
coefficients have been observed at 45° for bare tubes and at 
60° for finned tubes with respect to the horizontal bed axis. 
The heat transfer coefficient values obtained were in compari- 
son to the values obtained for tubes of the same configura- 
meomebut in a horizontal and vertical position [23]. 

3. Vertical Tubes 


The total heat transfer coefficient for a single ver- 
tical tube generally increases with a decrease in tube length 
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and diameter, a decrease in the fluidized particle size (sizes 
less than about 1 mm (0.039 in) at atmospheric pressure), and 
an increase in fluidization gas velocity up to u,- 
The heat transfer coefficient for a tube within a 
bundle increases with an increase in the tube spacing and 
is greater for tubes located nearer to the center of the 
bundle, at least for small bundles. 
The heat transfer coefficient for a finned tube 
generally increases with the fin height and a decrease in 
the frequency of fins along the tube length (in the range of 


21-10 fins/in). Heat transfer also improves with roughened 


and grooved tubes. 
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IV. EXPERIMENTAL APPARATUS AND PROCEDURES 


This section will describe the details and procedures as 
related to the experimental apparatus. The methodology and 


experimental procedures will also be discussed. 


A. DESCRIPTION OF THE FLUIDIZATION APPARATUS 

The fluidization container for the experimental investi- 
gation was constructed of one half inch thick clear plexiglas 
(Grade GM). The inside dimensions of the rectangular enclo- 
Sure were variable as one wall was capable of movement within 
the boundaries of the container. This movement allowed the 
Variation of the inside dimensions from six inches wide by 
Six inches deep by eighteen inches high to twenty four inches 
wide by six inches deep by eighteen inches high above the 
distributor plate. All joints were glued together and wood 
screws added in mounting the distributor plate for additional 
weight carrying Support. Mounted within each of the two side 
walls is a five inche wide by ten inch high strip heater 
that has the heating surface on the same plane as the inside 
plexiglas wall. Backing the strip heaters are two layers of 
insulation (silicone rubber sponge and plexiglas) which 
Brovide both rigidity and a mounting platform for the strip 
heater. Static pressure taps and bed thermocouple taps were 


installed at evenly spaced intervals as shown in Figure 27. 
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Pressure drops were measured uSing the static pressure taps 
with the output going to a Meriam Type W Model A324 sixty 
inch water manometer or a Model 33KA35 thirty inch water 
manometer via a pressure manifold. Air and bed temperatures 
within the experimental apparatus are measured with copper- 
constantan thermocouples as described in the EXPERIMENTAL 
HEAT TRANSFER PROCEDURE AND APPARATUS section of this thesis 
and recorded on a Newport Model 267B-TF2, Multipoint (12 
Channel) Digital Temperature Recorder. The complete equip- 
ment listing for this experimental study is given in Appendix 
=. 

The distributor plate utilized in the fluidization con- 
tainer was constructed of one half inch thick plexiglas 
(Grade GM) with a grid of 0.125 in diameter holes bored at 
0.25 in center to center spacing. To cause the required 
pressure drop which aids the air distribution and to prevent 
fae particles from falling into the air box, one layer of 140 
mesh stainless steel wire cloth was mounted on top of the 
distributor plate. 

The air waS Supplied by a Clements Cadillac Blower/ 
Suction Cleaner Model G12 with a nominal output of about 20 
scfm of air at a nominal line pressure of 1.5 psig. The flow 
Output data is based on experimental testing of the blower 
Since blower characteristic curves are not available. The 
air inlets to a two and one half inch inside diameter plexi- 


glas tube that has one half inch circular air ports spaced 
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@e One inch center to center intervals along its length. The 
meee ports are located at the underside of the tube allowing 
the air to exhaust downward into the air box. The plexiglas 
tube runs the length of the air box which has inside dimen- 
Sions of twenty four inches wide by six inches deep by six 
inches high. The air is unfiltered and prevented from exit- 
ing outside of the variable bed inlet area by an air baffle 
mounted to the movable wall above the distributor plate. The 
schematic line diagram of the air supply system is shown in 
Figure 28. 

The air flow rate was measured with a Fischer and Porter 
Company Model 10A3565A Flowmeter (rated at 11.1 scfm) which 
was corrected for local barometric pressure. Photographs 
showing two views of the overall experimental apparatus can 


be found in Figure 29 and 30. 


B. EXPERIMENTAL HEAT TRANSFER PROCEDURE AND APPARATUS 

Four different rectangular bed configurations using a 
Standard static bed height were investigated. For each bed 
configuration both an unheated and heated run were conducted 
with various pressure, temperature and visual data taken on 
Seen run. 

Once the bed size was adjusted, the bed of particles was 
initially fluidized at a point above minimum fluidization. 
When steady state conditions had been reached, the bed and 


air (inlet and outlet) temperatures were recorded anda 
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Figure 30. Photograph of experimental apparatus. 
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comment recorded describing the visual condition of the exter- 
nal boundaries of the bed. This was the starting point for 
the run. 

As a separate study, bed pressure drop was recorded as a 
function of several flow rates. This information was used to 
determine the point of minimum fluidization for the particular 
bed configuration and to serve as an indicator for the pre- 
sence of slugging or channeling conditions within the bed. 

After the base data run was recorded, power was applied 
to the two Watlow Silicone Rubber strip heaters by way of a 
Lamba Regulated Power Supply, Model LK345A. The power supplied 
to the heaters was computed by meaSuring the voltage and 
amperage. Energy losses by conduction were accounted for in 
the heater configuration and an energy balance computed for 
eaen Fun. 

Since the purpose of this experiment did not include a 
need to determine a heat transfer coefficient, provisions 
were not made to determine heater wall temperatures or the 
temperature gradients close to the heater wall. On all runs 
the expanded bed height covered the strip heaters ensuring 
that all heat transfer was either to the bed or losses 
througn the insulation. Radiation losses were neglected 
Since bed/wall temperatures were considerably less than 1000 
IK . 

At various heater power levels a series of temperature 


readings were recorded until flow characteristics were 
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determined for each bed configuration. The above experimental 
procedures were continued until all bed geometry configura- 
tions had been studied. 

The details of the heater configuration and mounting 
method are shown in Figure 31. Insulating end blocks were 
used for both mounting and end heat loss reduction for the 
two heaters. Heater insulation consisted of 0.5 in of sili- 
cone rubber sponge with a mounting block of 0.5 in of clear 
(Grade GM) plexiglas. Two insulated copper-constantan thermo- 
couples were applied to the outside of the plexiglas mounting 
blocks to measure their outer surface temperature. These 
temperatures were used to estimate the heat loss through the 
back of the heater by conduction for use in the heat balance 
of the model. 

The strip heater resistance element is constructed of a 
fine nicket alloy wire wrapped around a glass cord. The 
element wires are arranged in a parallel rectangular configur- 
ation thus confining any failure to a single element rather 
than the entire unit. Each heater is rated at 250 watts, 

115 volts with a temperature range of 0-200°F. A consider- 
ably higher temperature range is possible by cementing the 
heater to various metal plates, but is limited in this applica- 
tion due to the melting temperature of the plexiglas. 

To measure the bed temperatures a grid pattern of thirty 
retractable Omega Engineering ANSI T (Copper-constantan) 


twenty-four inch stainless steel sheathed exposed junction 
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thermocouple probes was used. By varying the depth of the 
individual thermocouples the boundary positions of the par- 
ticle and fluid flow cells could be determined and changes 


to the original pattern noted during the various runs. 


C. EXPERIMENTAL APPARATUS CALIBRATION 

Manometers used for static pressure drop and differen- 
tial pressure measurements were continuously rechecked for 
zero readings and readjusted as required. 

The temperature recorder and thermocouples were calibra- 
ted with a Rosemount Commutating Bridge Model 920A and 
Rosemount Calibration Bath Model 913A using liquid nitrogen 
and water. A least squares linear curve fit was applied to 
the resulting data for each thermocouple to obtain suitable 
constants for temperature correction. 

Technical data on the bed particles and experimental 
verification of the manufacturer's specifications can be 


found in Appendix F. 
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V. PRESENTATION AND DISCUSSION OF RESULTS 


Based on bed observations during the fluidizing process 
a very definite electrostatic charge built up within the bed 
container, see Figure 22. This phenomenon was characterized 
by the bed particles adhering to the walls of the bed after 
the air had been turned off and to the thermocouple probes 
during fluidization. Based on evidence provided by Miller 
and Logwinuk [49] this phenomenon will to some extent con- 
tribute to erratic results in the heat transfer and data 
measurements. Some of the factors wnich are thought to 
affect the magnitude of the electrostatic charge effects 
are (1) the dielectric constant of the bed material used, 
myeeche contact potential difference, (3) the ratio of fluid- 
ized bed particles to gas, (4) the properties of the gas, (5) 
the temperature of the solid-fluid system, and (6) the nature 
of the surface of tne bed particles. For the purpose of this 
thesis, the electrostatic charge existance will be acknow- 
ledged but no data correction will be considered warranted. 


Further study on this area is, however, considered worthwhile. 


A. FLOW PATTERNS 
The particle movement along the wall at any level is in- 
fluenced by the movement within the fluidized bed at that 


level, as it is the movement of the inner particles which 
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drag the outer particles along the wall. Therefore the wall 
velocity profile should represent approximately the average 
particle velocity at the various levels within the bed. Had 
the motion been smooth with all particles starting at the 

top of the bed and flowing straight down to the bottom of the 
bed then the ideal, Figure 2, flow pattern would have occurred 
and the velocity of the particles could be assumed to be con- 
stant at the wall with a positive linear gradient into the 
axis of symmetry for the bed. In actuality this was neither 
observed nor anticipated within the rectangular bed under 
study due to the geometry involved. 

hierieminltial Geometric configuration, that of a square 
@Belumn (6x6xl2.5 in), the particle motion appeared to conform 
memmeene particle circulation pattern of Figure 33. The fluid- 
ization somewhat resembled waves in that there were intermit- 
tent fronts moving smoothly up through the bed which terminated 
in a percolating (rupture) action at the top of the bed. The 
domed bulge, Figure 34, of particles was observed roughly in 
the center of the column and after rupturing, the bed particles 
would be thrown to the sides where they would start down the 
bed walls. 

It was not possible to observe any distinct points of 
emergence into the bed of the particles in the downward tran- 
Sit but a definite entry point was visible at the bottom of 
the bed. The characteristic shape of the non-movement area was 
elliptical along the bed bottom and curving up at the corners, 


See Figure 35. 
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Initially, the shape of the non-movement area was not 
Slliptical. If the side of the apparatus as shown in Figure 
27 is considered to be the front face then bubbles were ini- 
tializing in the front left corner of the bed along the 
movable wall. This snift in symmetry of the fluidizing cell 
generation was resulting in a shallower non-movement area at 
mae left front corner and a larger build up at the right 
Bront corner. The cause of the problem turned out to be 
Meee the distributor plate was not at a right angle to the 
air flow. Downward motion was observed at all four walls and 
the non-movement zone was located within one inch of the 
Bottom of the bed. 

It was readily observable that the velocity of the par- 
ticles emerging at the center of the bed was higher than that 
transiting down the face of the wall. This differential 
while not empirically measured during this study was attri- 
buted to frictional drag at the wall both mechanical and as 
a result of the electrostatic charge present at the bed wall. 
There 1S some justification provided by Toomey and Johnstone 
[62]. Marsheck [43] also gives an excellent treatise on the 
velocity profile across a gas-fluidized bed. 

The height of the domed bulge would pulse during the 
Fluidization more so as the air velocity was increased. If 
the air blower was fully energized with the initial bed con- 
figuration then the entire top of the bed would rise up as 


amass and then collapse on itself. This was an example of 
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slug flow where the rising bubbles had coalesced into one large 
bubble stretching from one wall to another. Slug flow was not 
achieved or expected in any of the other configurations as the 
bed was not deep enough nor was the air velocity high enough 

to provide the right conditions for the bed cross-sectional 
areas present. 

For the air flow that just permitted particle circulation 
the bed height fluctuations were negligible although the par- 
ticles at the wall still moved downward in an irregular 
Slightly surging movement. It is felt that at some point 
in the bed for this minimum fluidization velocity that the 
flow pattern internally is such that the flow pattern dis- 
tributes the particle velocities more evenly throughout the 
column. In other words, the circulating motion caused by 
the higher velocity in the center and slower velocity along 
the wall gives way to a more random or eddying motion that 
evenly distributes the small air bubbles within the two 
velocity boundaries thus preventing rapid bubble coalescence. 
These eddy cells once formed in the mixing region of the bed 
Will be moved up through the bed by the surges resulting 
from the cell rupture at the bed top. Between each surge 
line it is theorized that only one eddy cell will exist and 
that it will grow as it progresses through the bed until it 
coalesces with the gas in the central cell core at the bed 


EoD Prior to cell rupture. 
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It was also noted that as the thermocouple probes were 
lowered into the bed that a flattening of the bed top would 
occur. This effect became more pronounced as the probe pene- 
tration increased and is felt to be as a result of the probe 
Meoviding a grounding channel for bubble/particle flow. [In 
effect the probes acted as a calming agent in the agitation 
@e the bed without changing the overall bed flow cell 
configuration. 

When the movable wall was positioned for a 6.25 inch width 
me = 37.5 in*) and also for an 8 inch width By Sas mo) the 
bed expansion observed at the minimum fluidizing air velocity 
was 0.125-0.375 in. When the bed width was then increased 
to 10 inches (A. = 60 in*) and subsequently 12 inches (AQ = 
72 in*) the bed expansion observed did not exceed 0.25 in. 
This measurement does not take into account the height in- 
Crease of the cell dome prior to rupture. There was no 
general horizontal movement of the particles near the wall 
except at the top of the fluidized bed and in a rather 
Burbulent mixing region near the bottom of the fluidizing 
Bed. 

As the movable wall was repositioned thus increasing the 
bed air inlet area the bed became an aggregative fluidized 
bed. This resulted in multiple cell formations within the 
bed vice the single cell observed in the initial bed config- 
mration. Unfortunately, there was no observable pattern 


to the cell formations as their overall positioning was 
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random in nature, see Figures 36 and 37. Only one cell was 
observed to break the surface of the bed at any given instant. 
Occasionally, small bubbles were noticed breaking the surface 
of the bed between the larger bubble cells. This was thought 
to be caused by channeling between the larger bubble cell 
boundaries. 

There was no noticeable change in the outward appearance 
of the particle cell flow patterns in any of the bed configura- 
tions when the external heaters were energized during 


mPatdization. 


See FE RESSURE DISTRIBUTION 

Static pressure readings were taken at various intervals 
along the height of the fluidized bed on all experimental runs. 
In order to determine the gas mass flow rates for minimum 
fluidization, measurements of the pressure drop across the 
bed were also made at various flow rates. All data obtained 
during the various runs can be found in Appendix H. The mass 
Flow rate at the maximum pressure drop was then taken as the 
point of minimum fluidization in each case. The maximum 
pressure drop was used because it is the most sharply defined 
point on the plots, see Figure 38. 

The results of plotting the pressure drop versus gas mass 
velocity profiles produced very similar patterns which are in 
general agreement with the expected curve for the distributor 
design used. As can be seen from Figure 38, when the air inlet 


bed cross-sectional area was increased the curve would shift 
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Mmagume 30. Bed pressure drop vs. gas mass velocity. 
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to the left. This displacement may nave been as a result of 
the lessening in the contribution of the mechanical wall 
friction drag to the total overall drag experienced by the 
bed with the increased air flow for fluidization. Neverthe- 
less, it would appear as the bed air inlet cross-sectional 
area was increased the gas mass velocity required for fluid- 
ization would also increase. Additionally, in all cases ob- 
served, the log Ap value for minimum fluidization was 
approximately 1.25. 

The results of the pressure measurements seemed to indicate 
that mechanically smooth fluidization was taking place. 
Channeling and/or slugging was not occurring at the minimum 
fluidization operating conditions although it was possible 
to demonstrate both channeling and slugging for the initial 
bed configuration. It was noted that during fluidization 
small fluctuations were visible in the liquid levels in the 
manometer that was measuring the static pressure in the bed. 
This operating characteristic is in agreement with similar 
Conditions observed by Parent et al. [52] in their study of 
gas fluidization. 

The chain of events in measuring the pressure drop is as 
follows. Once the bed was initially settled the air flow 
Was initialized to the bed. As the air flowed through the 
bed, a corresponding pressure drop occurred. This pressure 
drop increased linearly with the gas flow until at some point 


the bed began to expand. At this point the pressure drop 
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became nonlinear yet still increased to some maximum value. 
This maximum value was picked as the minimum fluidization 
point as previously stated. The pressure drop would then 
decrease until the bed started to circulate or fluidize at 
which time the pressure drop would stabilize at some point 
lower than the maximum value recorded earlier. The pressure 
drop remained steady as the air velocity was increased further 
but soon started to rise again. This new rise inthe pressure 
drop across the bed was an indication that slugging was start- 
ing to set in and was especially prevalent in the initial bed 
geometry. It was also noted that this secondary pressure 

drop increase was less prevalent as the bed cross-sectional 


area was increased. 


Pet aeMPERATURE DISTRIBUTION 

Air entered the air box at temperatures between 70° and 
120°F and in every case left the bed at a temperature within 
1°F of that of the top of the bed, as measured at a point in 
the axial center of the bed. Only two data runs were made 
with the bed instrumented to record the temperature distribu- 
tion. The first run (A, = 37.5 in*) had an average air inlet 
temperature of 108°F +2.5°F and attained a constant air outlet 
temperature of 125°F +0.4°F, see Figure 39. The second tem- 
perature instrumented data run (A. =i naa aneaverage 
air inlet temperature of 116°F +2°F and attained a constant 
air outlet temperature of 129.5°F +0.7°F, see Figures 40A and 
40B. 
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In the graphics of Figures 39 through 40B it was necessary 
to plot the profiles for each thermocouple probe as there were 
too many variables to produce a true three dimensional temper- 
ature profile. The x-axis of each graph is the height of the 
probe tip in inches above the bed distributor plate. The y- 
axis of each graph is the temperature in °F. A probe grid is 
provide for each bed configuration with the probe number for 
reference. In all cases the temperature recorded by each 
probe above the 1.5 inch height was within the constant air 
Outlet temperature range specified above for each data run. 
All temperature data recorded can be found in Appendix H. 

As can be readily observed, the bed quickly attains a 
constant temperature except for a very small region at the 
bottom of the bed and a narrow layer near the container walls. 
The area near the walls is by assumption as the experimental 
apparatus was not instrumented to read that area. 

It appears that the particles were heated by the heater 
as they transited TINS kecheehe Surface. Then, when they reached 
the bottom of the bed, thermal equilibrium with the incoming 
cold air was attained almost instantaneously. Except for the 
layer along the walls of the bed and in the small turbulent 
mixing region at the bottom of the bed, thermal equilibrium 
between the particles and the gas very nearly exists through- 
Out the fluidized bed. 

Since the inside walls of the bed container were not 


instrumented with thermocouples it was not possible to 
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determine if there existed a temperature differential between 
the walls which did not have heaters and the bed bulk 
memperature. 

In order to determine the effect of the air velocity on 
the temperature readings of the thermocouple probes several 
runs where conducted at different air inlet velocities with 
the probe tip placed six inches above the top of the distri- 
butor plate. Small temperature fluctuations of at most 1°F 
were observed while the bed was fluidized and being heated 
under steady state conditions. The air flow was then ter- 
minated while monitoring the bed temperature readings. It 
was observed that the thermocouple reading would initially 
Stabilize at a temperature within the range of the previous 
fluctuations and did not change for several seconds. A 
gradual rise then would commence due to the heat conduction 
from the heaters in the walls. Since the temperature rise 
Mid not exceed 1°F at the end of four minutes it would seem 
to indicate that the bed thermocouple readings were basically 
independent of the fluid velocity. It would also tend to 
Support the assumption made earlier in this thesis that the 
particle temperature and gas temperature at any point in the 
bed can be considered to be at equilibrium. 

The only change in the mixing cell pattern within the 
fluidizing particle cell appeared to be the production of 


the temperature gradient at the bottom of the bed and an 
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increased bulk bed temperature. This was in comparison to 
a steady state bulk bed temperature equal to the air inlet 


temperature prior to energizing the plate heaters. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


= CONCLUSIONS 

The following conclusions may be made from the results 
of this investigation: 

1. Upon fluidization of the bed an electrostatic charge 
will build up within the bed. Although no data correlation 
error is used in this thesis, there most probably is some 
effect on data results. 

2. With a square column bed the particle flow cell formed 
is basically as shown in Figure 33. Only one major cell is 
Bermed with its axis of ascent roughly located at the center 
meechne bed. 

3. With an expanded rectangular bed several particle 
flow cells form in a random fashion. There is no pattern 
meme Overall cell coverage and the bed can be termed a 
bubbling bed. 

4. That since the bulk bed temperature remains constant 
throughout the bed that there must exist eddy cells within 
the bed that continually mix the bed particles as they move 
up through the bed. Otherwise, a gradient would form between 
the heater walls and the center of the bed. 

>. That the velocity in the center of any cell was 
greater than the velocity of the particles being forced down 
the walls. In the expanded bed the particles moving down the 
Walls were a mixture of the various ruptured cells. 
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6. That mechanically smooth fluidization was not taking 
place at all times. Occasional instances of channeling were 
present and could be distinguished by small pertubations 
to the bed surface that did not change position - a simmering 
appearance. 

7. That the distributor plate used did provide enough 
of a bed pressure drop to allow fluidization to occur but that 
the tendency for channeling or slug flow was high if the air 
velocity was increased much past the minimum fluidization 
velocity. 

8. That a temperature equilibrium between the incoming 
fluidizing air and the bed particles was reached almost in- 
Stantaneously in a small mixing region of the bed located at 
the bottom of the column. 

9. That the only temperature gradient that existed in 
the bed was located at the bottom of the bed where the mixing 
region was located. The remainder of the fluidized bed stayed 
at a constant bulk bed temperature irregardless of the number 
of flow cells present. 

10. That the air velocity did not affect the bulk bed 
temperature readings. It does, however, affect the depth of 
the mixing region at the bottom of the bed. 

ll. That energizing the plate heaters did not visibly 
alter the particle cell flow pattern but that it did produce 


a gradient in the temperature profile of the bed interior. 
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fi2. That the rectangular bed configuration under study 
would trip into fluidization when the log Ap across the bed 


epproached 1.25. 


B. RECOMMENDATIONS FOR FURTHER STUDIES 

Based on the results of the experimentation to date, 
research conducted in preparing this thesis and the problems 
experienced in trying to get the fluidization apparatus to 
work the following areas of study are recommended for future 
thesis study: 

1. A determination of what effect the electrostatic 
charge which builds up during the fluidization process has 
on the heat transfer or flow cell formation. 

2. A study of the resulting changes in the cell formation 
patterns when bed particles of sizes larger than 1 mm are used. 
Particles in this size range are of much more practical use in 
the metallurgical industry. 

3. A study of the flow cell patterns with beads smaller 
than 300 microns which should give enhanced mechanical 
mruidization. 

4. A study of the effect of particle stratification 
(elutriation) when particles of different sizes are combined 
in the bed. The combination can be of varying sizes/rough- 
ness/sphericity or any other physical property of the bed 
particles which would promote the separation of the particles 


Mnto distinct layers during fluidization. 
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5. A review of the cell pattern formation when using 
particles from various manufacturers (same or different bead 
materials) while maintaining a constant particle size. 

6. Conduct a pressure study on the bed with different 
distributor configurations. This would include such designs 
as porous plates, baffles, calming chambers, filter screens 
Or various combinations of the aforementioned designs. 

7. Study the effect of changing the location of the heat 
source from opposite wall external heaters to other configura- 
tional combinations such as all four walls, adjacent walls, 
Single wall, submerged heater, etc. 

pepeeconauct a velocity profile study similar to that con- 
ducted in Reference [43]. That study used a modified thermis- 
tor anemometer probe to determine the flow pattern within a 
eeecular column bed. 

Additionally, prior to conducting future research the 
following apparatus modifications are suggested: 

l. Replace the left side sponge rubber gasket on the 
movable wall with a new gasket. 

2. Implant a pair of thermocouples in the unheated walls 
of the fluidized bed. 

3. Mount a pair of thermocouples on the surface of each 
strip heater to enable the monitoring of the heater wall 
temperatures. 

4. Manufacture a stand for the cabling to the thermo- 


couple probes. The stand should be of sufficient height 
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mo prevent the bending of a probe from the weight of its 
cable when the probe is retracted from the bed. 

5. Change the mounting method of thermocouples 35 
through 38 so that they are mounted in the backing plate 
vice on the backing plate. 

6. Unless the area of study is to focus on the electro- 
Static charge built up during fluidization, provide some 
means to ground out this phenomenon. 

7. Provide a leveling system for the fluidization appara- 
tus that will ensure that the fluid flow direction out of 
the distributor plate is parallel to the bed container walls. 

8. The probe taps on top of the fluidizing apparatus 
snould be bored to 5/64 inch for ease in inserting the probes 
during data runs. Additional nylon bushings should also be 
ordered. 

9. A second flowmeter of greater capacity than 11.1 scfm 
Should be installed in the air supply line. This will be re- 
quired to measure the higher air velocities needed to fluidize 
the bed with further expansion. 

10. A cooler air source should be used in future runs to 
enhance the temperature differential at the bottom of the bed 
during fluidization with the wall heaters energized. 

ll. A new cap for the fill vent is requird to prevent the 
glass beads from getting between the cap and plexiglas thus 


jamming the cap. 
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C. CONCLUDING REMARKS 

Particle motion is a complex function of the particle and 
Bes properties, fluidization velocity, gas distributor and 
geometry of the bed. While a simple correlation can easily 
be made it must be tempered in its general application by the 
fact that it will most likely not give adequate allowances 
for the way in which the particles circulate and subsequently 
change with changes in the bed particle material, fluidizing 
conditions, scale and design of the fluidization apparatus. 

The interpretations presented in this thesis are based 
on the premise of visible particle movement at a given bed 
height and must be qualified by the fact that they represent 
a very limited amount of data. 

There is a general Opinion among researchers that the 
bubbles of air pasSing through a fluidized bed grow in size 
by coalescence of smaller bubbles into larger bubbles. Even- 
tually given a deep enough bed the bubbles should coalesce 
into one large bubble located in the center of the bed in 
order to maintain the uniformity of the bed [43]. While 
this Opinion has been generated as a result of experimenta- 
tion uSing circular beds it would appear that it can be 
equally applied to rectangular beds. No doubt had the 
expanded rectangular bed been deep enough, eventually it also 
would have conformed to the aforementioned opinion. 

The determination of particle flow patterns in a fluidized 


bed would be of extreme value in the establishment of a kinetic 
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or mass transfer model of a fluidized bed. The fact that 
there are a multitude of interrelated variables involved, 
unfortunately, has made the modeling a theoretical tool vice 
the needed operational scale-up tool required for industry. 
More study will undoubtedly be required until a workable 
industrial model can be produced. 

It has also been an observation that while the majority 
of the fluidization research was conducted by American 
researchers in the the late forties and early fifties, that 
the bulk of the recent research has been conducted by the 
Russians with some notable exceptions. This trend is neither 
understood nor considered warranted by the researcher con- 
Sidering the very pronounce heat transfer improvement demon- 


strated through the use of fluidization. 
& 
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APPENDIX A: NONDIMENSIONAL CORRELATIONS 
FOR EXTERNAL WALLS 
These correlations are but a sampling of the many 
fluidization correlations for external walls that have 
been theorized. 


1. Investigator - Baerg et al. [3] 
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Investigator - Das and Sarkar [11] 
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Investigator - Dow and Jakob [13] 
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Investigator - Gamson [17] 
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Investigator - Lemlich and Caldas [35] 
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Investigator - Leva [36] 
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Investigator - Levenspiel and Walton [41] 
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meee investigator — Mickley and Trilling [48] 
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12. Investigator - Rao and Kaparthi [56} 
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13. Investigator - Richardson and Mitson [57] 
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Remarks - Liquid-fluidized system with 
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14. Investigator - Rukenshtein [58] 
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15. Investigator - Toomey and Johnstone [62] 
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APPENDIX B: NONDIMENSIONAL CORRELATIONS 
BOn seve Roe SODIES 


These correlations are but a sampling of the many 


fluidization correlations for immersed bodies that have 


been theorized. 
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Investigator - Andeen and Glickman [2] 
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Investigator - Gelperin and Einstein [18] 
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Remarks - Vertical tube bundle. 
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Investigator - Gelperin et al. [20] 
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Investigator - Miller and Logwinuk [49] 


Coprelation — 


Ges k 050.7 -2 k 2 4 
S £ 





he — las i - 
d 0. CO - O36 
p oo = 
Remarks - Vertical tube bundle. 
Investigator - Petrie et al. [54] 
Correlation - 
jeeeG 7 
hi Den Gc | 0.33 Dr 
| Eee | 
Ke Sie qd, 
Remarks - Staggered and in-line tube bundles. 
Investigator - Sarkits et al. [60] 
Correlation - 
0.45 0.16 
Cs D 
Mie — On Olesmre* Ga-’ pros? | =— — 
Ce qd. 
Remarks - Coil. 
Investigator - Traber et al. [63] 
Correlation - 
0.1 3 0.16 
D a 
i momo O4mee- 2 ea opr? | _— 
da d 
Pp Pp 


Remarks - Coil. 


eZ 





JLo 


nS 


Investigator - Vreedenberg [69] 
Correlation - 


@ 0.3 5 


Nu = a 
arte Vin 


Remarks - Tube position varies. 


Investigator - Vreedenberg [70] 


Correlation - 


ee NaS 


Remarks - Horizontal tube with- | ————| < 2050. 


Correlation - 


Gabe 0.3 0 


Remarks - When |————| is greater than 2050 but 


U Pr 


less than 2550 the value of Nu will be 
the average of the two Nu correlations 
above. 
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14. Investigator - Wender and Cooper [74] 
Correlation - 


0.4 3 0.8 


Nu = 0.0119 (1-+ ¢) )_—<c a = 


Remarks - Vertical tube. 


Pree livescigatos — Zenz and Othmer [80] 


Correlation - 
Woe) 6:7 
G Din 
Nu = 14 |—— Sree tet || oreall 
Caf Be 
Remarks - Horizontal tube. 





meee Noik Cs NONDEIMENSIONAL CORRELATIONS FOR THE 
MAXIMUM HEAT TRANGPER COEFFICIENT 
These correlations are but a sampling of the many 
fluidization correlations for maximum heat transfer coef- 
ficients that have been theorized. 
Meee ifivestigator - Baerg et al. [3] 
Correlation - 


Or O00 27 Pop 


homax = 49 log ; 


p 


Remarks - None. 


NO 


Investigator - Chechetkin [8] 


Correlation - 
0.8 0.4 
OC On 7 RE oot Px 
Nu = 
max gq 069 
D 
Remarks - Re =m Woe a > - . 
Se opt 


Dee investigator - Chekansky et al. [9] 
Correlation - 


0.04 D ~ 0.1 2 


A 
h = 28.2 0 Uae k 0.6 a 0.36 
wmax : S fe p 


Os 
O 


Remarks - Staggered tube bundle. 


Ws 





Investigator - Gelperin et al. [19] 


Correlation - 
0.2 5 
Dr Dip 
Nu Om teGaee | 1 — ——| | + 
oon Zit 
Y T 
Remarks - Staggered horizontal tube bundle. 


Investigator - Gelperin et al. [21] 


Correlation - 


Remarks - Vertical tube bundle. 


Investigator - Gelperin et al. [22] 


Correlation - 
0,2 5 
Din 
Nu = jewomea- 4 )| 1 = == 
ax 
vi 
Remarks - In-line horizontal tube bundle. 


Investigator - Jacob and Osberg [30] 


Correlation - 
a= ae (awe Oe 1 =nexp» { Spe5! 
Remarks - AY and p determined by experimentation. 





JUDG 


ite 


Investigator - Maskayev and Baskakov [44] 


Gorrelatien). — 
Nu = (0), 21, @erte? 
m 


Remarks - Vertical cylinder in coarse particle bed 
ieee alm Cae <sosclO - 


Investigator - Pillai [55] 


Correlation - 


Nu =——OmsG5T> ° - Gals 


Remarks - Sphere with Bi < 0.04 
Investigator - Ruckenstein [58] 


Correlation - 


Nu = Re 0.42 3 Gat Pr ?3 3 
max Stic 


er = 0.5 8 
Remarks eee 0.09 Ga : 


Investigator - Sarkits [60] 





Correlation - 
045 O21 6 0.45 
ps D ae 
me 0.4 2 0.33 in She 
Nt aaa = 0.0087 Ga Pr 
S d d 
pf Pp p 
Remarks - Re opt Ome Ga | lamigar cegion) - 
Correlation - 
0.1 G13 0.16 
“ps D are 
Nu = (5 (0 LS) (ey ite le — — 
max 
S d d 
pt Dp p 
Remarks - Re ont =—OMuiGmGaten (EM rbu lent region) . 


130 





2s 


ise 


a. 


Investigator - Traber et al. [63] 


Correlation - 


Nu Onde! Ga. Pree? | —— eaae 


ss 

my 

% 
Qu 
Qu 


Remarks - Re. =m 55 Ga. 


pt 


Investigator - Varygin and Martyushin [67] 


Correlation - 


Nu ="(eeG Ga 
m 


Remarks - Re =i | eee ?:> : 


pt 


Investigator - Zabrodsky [78] 


Correlation - 


= Ole 0.6 — 0.36 
Danese = O2od Bij Naas Ce 


Remarks - None. 
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APPENDIX D: GENERAL EQUATIONS AND UNITS 


This appendix will be broken down into two sections. 
The first section will cover the general nondimensional heat 
transfer/fluids numbers used in this thesis. The second 
section will annotate the units to be used with each equation 
that is numbered within the main text. This requirement was 
necessitated by the use of different units (same variable 
nomenclature) among the various researchers cited in the 
thesis text. 


A. GENERAL NONDIMENSIONAL HEAT TRANSFER/FLUIDS NUMBERS 


1. Reynolds Number (Re) 


PeeeeGalileo Number (Ga) 


3 = 
qd ° Pe (2 Pe VY 
Gas: —___.. 


rS2 


Nusselt Number (Nu) 





Stanton Number (St) 


W 
Sst = 
u 


Sie Pe Yo 


Biot Number (Bi) 


B. EQUATION UNITS 


Ie 


Equation (3.1) 


1, = Way ats 
1 -~ ft 
e 
=a - dimensionless 
Die, = Weiiy,/sere 
Oe ~ lbm/ft? 


Heuation (5.2) 


Ent ~ dimensionless 
V - ft? 
Pp 
ie lb/ft? 
L af - ft 
A - ft* 
Cc 
0 - lbm/ft? 


Ls 





P¢ - lbm/ft? 


Bomacton (3.3) 


= = 2 

Gf lem nc—£t - ) 
D =e 

p 
u - centipoise 
= lbm/f£t? 
an lbm/f£t? 
Equation (3.4) 
D = fie 

p 
D. - in 

1 
x - dimensionless 
Equation (3.5) 
Re wf - dimensionless 

a ee 

G Af lbm/ (hr-ft*) 
D - in 

Dp 
u - centipoise 


Equation (3.6) 


a at ie? 
Gat lbm/ (hr-ft*) 


meations (3.7) 


n - dimensionless 
Ge - lbm/ (hr-ft?) 
a) | lbm/ (hr-f£t7) 


Bavation (3.8) 


1S - dimensionless 
m~ - dimensionless 
G- - lbm/ (hr-ft?) 





HO’ . 


a. 


IL vas 


G ~ lbm/ (hr-ft?7) 


mf 


BawWat lon. (3.9) 


AR ~ 


x 
! 


12) 
! 


Q 
| 


e - 


h = 


1G _ 


Dow 


wn 
{ 


Q 
j 


H) 
t 


EMBO, qiia 
BTU/ (hr-ft*- °F) 


ee 
Ln 


Bewation (3210) 


~ BTU/ (hr-ft7-°F) 


BTU/ (hr-f£t- °F) 
lbm/ft? 
BTU/ (lbm- °F) 


Sec 


Beuation (3.11) 


BTU/ (hr-ft?- °F) 
BTU/ (hr-ft-°F) 
lbm/f£t° 

BTU/ (1lbm- °F) 
Be7 ar 


fete 


Bauation (3.12) 


BTU/ (hr-ft7- °F) 
BIU/ (hr—f£t- FP) 
BTU/ (lbm-°F) 


l/fsec 


LoS 





Device 


Temperature 
Recorder 


Temperature 
Recorder 


Thermocouple 


Water 
Manometer 
0 in) 
Water 
Manometer 
60 in) 


Flowmeter 


Variac 

Air Blower 
Regulated 
Power Supply 


Thermocouple 
ealLibration 


SErip 
Heaters 


Weight Scale 


Spotical 
Microscope 


APPENDIX E; 


Manufacturer 


Newport 


Newport 


Omega 
Engineering 


Meriam 
Instrument 


Meriam 
Instrument 


Fischer and 
Porter 


General Radio 


Clements 
National 


Lambda 
Electronics 


Rosemount 
Engineering 


Watlow 
Electric 


Mettler 
Instrument 


Bausch & Lomb 
Instrument 


it So) 


BROUReEMENT  LEStTiINnG 


Model # 


26 7B-TF2 


267B-TF2 


ANSI T 
Type S 
Sheathed 


Type W 
33KA35 


Type W 
A324 


LOA3565A 


Type V10 


GiZ 

fies .4 on M 
Bridge-920A 
Bath-913A 
250 Watts 


Type H15 


31-32-45 


Serial #4 


S200 39 2 


So200293 


MG Jos 


44F554 


6908A2165AlL 


Tao S6 


Crete 2a: 


JEG, 


8044 


ESI 23 


B=2o 07 ks 





Pecos PARTICLE ANALYSIS 


The particle specifications provided by the manufacturer, 
Potters Industries, are as follows: 
l. Physical Properties 
a. Size - Minimum 90% by weight in range on smaller 
U.S. Sieve, maximum of 10% by weight on larger U.S. Sieve. 
Tested in accordance with ASTM-D-1214-58 and MIL Spec. 
G-9954A. 


b. Free Iron Content - Maximum of 0.1% by weight. 


c. Chemical Content - Soda-lime silica glass. The 
solid glass spheres contain no free silica. 


Greco becifve Gravity = 2.45-2.55. 
e. Hardness - DPH 100g load-515Kg/mm?. 


f. Color - Clear and colorless, free from surface 
ems: . 


q. Roundness - Minimum percent round by roundometer 
is 90%. Determined in accordance with ASTM-D-1155-53 or 
foeeopec. G-9954A. 


h. Broken or Angular Particles - Maximum of 3% by 
count. Tested in accordance with MIL Spec. G-9954A, 


1. Air Inclusions - Not more than 10% of the beads 
Shall show air inclusions of more than 25% of their surface. 
Tested microscopically immersed in 1.5 refractive index 
falta. . 

4. Average Settled Bed Density - 96 lbm/ft?. 

k. Average Particle Density - 156 lbm/ft’. 

imecocrtriecrent Of Friction (Static) -10.9-1.0. 

mee thermal Properties @ 70°F 


a. Specific Heat - 0.180 BTU/(lbm-’F). 


ey 





De LiommaleConductivity = 0.515 BIU/ (hr=tti—°F) . 
3. Belectkrical Properties @ 77°F 

Sere lecteric Gonstant (1 Kc) - 7.6. 

Dee DeewvolimemResistivyaty ~ 6.5X10>* Ohm-cm. 
4. Size Analysis (Potters P-020 Designation) 

a. Basic Sieve Analysis - U.S. Sieve 35 to 50. 


byeeeoverage Particle Diameter - D_ = 0.0122 inch 
LO microns). P 


c. Maximum Percentage True Spheres - 90%. 
Photographs of two groups of the particles are shown in 
Figures 41 and 42. Upon microscopic examination the beads 


appeared extremely spherical in nature with some air 


inclusions present in all beads. The inclusions appear as 
black areas on the photographs. There was little variation 
in particle diameters between the two sample groups. The 


diameter for each group was taken as the arithmetic mean of 
Microscopic meaSurements on a representative sample of 50-60 
beads in each group. Group I had an average particle 
diameter of 0.0111 in while Group II had an average particle 
Giameter of 0.0110 in. For the purpose of calculations in 
this thesis an average particle diameter (D,,) Seeoweell in 
was used. 

A basic sieve analysis was performed on two samples of 
known volume and weight to verify the manufacturer's 
Specifications. The sieves where U.S. Tyler Company stand- 


ard screen mesh of 48 and 60 mesh size. Less than five 
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meoure 41. Micropnotograon of beads-Group I. 














gy. 
She Recipe ty 
4 





Senate 


Figure 42. Microphotograph of beads-Group II. 


29 





percent by weight of the sample particles was larger than 
the 48 mesh screen and less than one percent by weight of 
the sample was smaller than 60 mesh. The Tyler screen scale 
of 48 (U.S. Sieve No. 50) has a sieve opening of approximate- 
ly 0.0117 in and the Tyler screen scale of 60 (U.S. Sieve 
No. 60) has a sieve opening of approximately 0.0098 in. 
Water displacement of a known weight of beads was used 
to determine the average particle density. The value 
Seeained Was 154.6 Ibm/ft* which is considered to be within 
acceptable agreement of the manufacturer's specifications. 
Likewise, the average settled bed density was obtained by 
measuring the volume occupied by a packed bed of beads of 
known weight. The value obtained, 93.2 lbm/ft*?, is also 


within acceptable limits. 
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APPENDIX G: EXPERIMENTAL UNCERTAINTY ANALYSIS 


The source of error in the results presented is consid- 
ered to be due to instrument precision and accuracy, 
inaccuracies in geometrical measurements and inaccuracies 
in data recording. The following is a breakdown of the 
experimental uncertainties. 

1. Distance Measurements 

The tape attached to the side of the fluidizing 
PPparatus has a minimum scale division of 0.25 in. Distance 
readings were made with this scale to the nearest 0.0625 in 
for both the bed height and width. While it was possible to 
jet static readings with a precision of +0.5%, the expanded 
bed heights were more on the order of +2% with a worst case 
precision of +4%. 

The depth measurement of the thermocouple probes was 
based on measuring the removal distance of the probe from 
its tap at the top of the apparatus. The smallest scale 
division for this withdrawal measurement was 0.0625 in. 
However, since the probes were not rigid structures, the 
bending in the probe lengths resulted in precision estimates 


=i +2 


fe)\ 


2. Particle Diameter Measurements 
The optical microscope used to measure bed particle 


diameters in the two samplings of beads has a minimum scale 
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division of 0.0001 in which allowed reliable measurements up 
to 0.00005 in at 100X magnification. The diameter precision 
based on scale readings is +0.5%. However, operator error 
in reading the diameters places the precision estimate more 
Bealistically at +1%. The standard deviation in particle 
diameter for the two bead groups was 0.00042 in and 0.00045 
on). 
3. Particle Density Measurements 
The scale used to measure the particle bead weights 
has a minimum division of 1 mg which equates to a precision 
of +0.6%. The minimum volume scale was 1 ml which equates 
to a precision of +1%. Due to spillage in filling the 
measuring containers with water and beads in addition to the 
losses during the experimental runs it 1S estimated that the 
error in density measurements is closer to +3% overall. 
4. Temperature Measurements 
Although the temperature recorder gives readouts in 
tenths of degrees, after taking into consideration the 
fluctuations experienced during fluidization and round off 
errors the overall error in temperature measurements 1S 
Pomsidered to be +1° F. 
9. Pressure Measurements 
The minimum manometer scale division was 0.1 in H:20 
with readings of 0.05 in H20 possible for a full scale 
precision of +0.2%. While this precision was possible during 


static bed conditions it was not possible to be as accurate 
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once fluidization occurred as the pressure level would 
fluctuate within the manometer column. For the worst case 
of the pressure measurement in the column, the pressure 
fluctuated +0.3 in H20 over the maximum pressure recorded 
Bor a recorded error accuracy of +£2%. 
6. Flow Measurements 
The minimum flowmeter division was 1% of 11.1 scfm 
Mee. llil scfm. This gives a scale precision of +1.2% for 
the range measured. Precision in reading the float position 
especially with fluctuations during fluidization caused 
additional inaccuracies. As a result, the overall flow 
measurement was probably within 15%. 
7. Miscellaneous Measurements/Calculations 
The following uncertainties were determined for 
various calculations/measurements required for data reduction 


by the method of Kline and McClintock. 
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APPENDIX H: EXPERIMENTAL DATA 


The following temperature and pressure data was recorded 
during experimental runs of the fluidized bed apparatus. 
1. Pressure Data - Run #1 
a. Static Bed Height - 12.5 inches. 
Be otacrc Bead Width —- 6.25 dnches. 
Eee Veena Gress -Secitional Area = 37.5 in’. 
d. Bed Pressure/Gas Mass Flow Data - See Table l. 
e. Comments - 
(1) Channeling observed to start at the edges 
Of the bed at an air flow of 8.88 scfm. 
(2) Fluidization of the bed commenced at an air 
mbew Or 9.66 scfm. 
(3) Bed expansion during fluidization was 
oeez5—-0.375 inch. 
2. Pressure Data - Run #2 
a. Static Bed Height - 12.5 inches. 
be otatic Bed Width - 8.0 inches. 
Poem rnilet Cross-Sectional Area = 48.0 in’. 
dad. Bed Pressure/Gas Mass Flow Data - See Table 2. 
e. Comments - 
(1) Channeling observed to start at the edges 
of the bed at an air flow of 9.44 scfm. 


(2) Fluidization of the bed commenced at an air 
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mow Of 11.32 scfm. 


wer Be@eexpansion during fluidization was 


aoe 5—0.375 inch. 


3. Pressure Data - Run #3 
demo eatiec Bed Height —- 12.5 inches. 
b. Static Bed Width - 10.0 inches. 
Cn ano tN Gross -Sectelonal Area — 60.0 in*. 
d. Bed Pressure/Gas Mass Flow Data - See Table 3. 
e. Comments - 


(1) Channeling observed at the corners of the 


meat an air flow of 10.55 scfm. 


(2) Channeling observed along the walls and 


occasionally in the bed at an air flow of 11.66 scfm. 


(3) Fluidization of the bed commenced at an air 


mitew Gf 12.7/ sctim. 


(4) Bed expansion during fluidization was 0.25 


meh . 
4. Pressure Data - Run #4 
weeesteatde. Bead Height — 12.5 inches. 
b. Static Bed Width - 12.0 inches. 
Pe oer tEGross=oectional Mrea - 72.0 in’. 
a. Bed Pressure/Gas Mass Flow Data - See Table 4. 
e. Comments - 


(1) Channeling observed at the corners of the 


bed at an air flow of 11.66 scfm. 


(2) Channeling commenced in the main part of 
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the bed at an air flow of 12.77 scfm. 
(3) Fluidization of the bed commenced at an air 


mow Of 15.00 scfm. 


(4) Bed expansion during fluidization was 0.25 
mach . 
5. Temperature Data - Run #1 
a. Static Bed Height - 12.5 inches. 
b. Static Bed Width - 6.25 inches. 
See eon mriCross Sectional Area — 37.5 in’. 
We Veragc e aseeinlete temperature — 109° F +2.5° F. 
Seeeconctant) sin Outlet Air Temperature = 125° F 
m4” FF. 
f. Temperature Distribution Data - See Table 5. 
6. Temperature Data - Run #2 
a. Static Bed Height - 12.5 inches. 
b. Static Bed Width - 12.0 inches. 
eee eee tCross-Sectional Area -§/2.0 in’. 
d. Average Air Inlet Temperature - 116° F +2.0° F. 
SwCGOotant 71 Ourlet Air Temperature - 129.5° F 
Boe? F. 


f. Temperature Distribution Data - See Table 6. 
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